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GENERATION AND FLUID DYNAMICS OF 
SCATTERING AEROSOLS IN LASER DOPPLER VELOCIMETRY 

M. K. Mazumder, B. D. Hoyle and K. J. Kirsch 


SUMMARY 

An experimental study on the particle-fluid interactions of scattering 
aerosols was performed using monodisperse aerosols of different particle 
sizes for the application of Laser Doppler Velocimeters in subsonic turbulence 
measurements. Particle response was measured by subjecting the particles to 
an acoustically excited oscillatory fluid velocity field and by comparing the 
ratio of particle velocity amplitude to the fluid velocity amplitude as a 
function of particle size and the frequency of oscillation. Particle velocity 
was measured by using a Differential Laser Doppler Velocimeter. An aspirator 
type atmoizer was used to generate aerosols from a dilute suspension of 
monodisperse polystyrene latex particles. A Rapaport-Weinstock generator 
was used to generate aerosols of DOP droplets. The test aerosols were fairly 
monodisperse with a mean diameter that could be controlled over the size 
range from 0.1 to 1.0 micron. Experimental results on the generation of a 
fairly monodisperse aerosol of solid particles and liquid droplets and on the 
aerosol response in the frequency range 100 Hz to 100 kHz are presented. The 
present study indicates that a unit density spherical scatterer of 0.3 micron- 
diameter would be an optimum choice as tracer particles for subsonic air 
turbulence measurements. 
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I. INTRODUCTION 

In a transonic or supersonic flow field where the flow medium is seeded 

with a scattering aerosol for velocity measurement with a Laser Doppler 

Velocimeter, a considerable amount of velocity lag may exist between the 

fluid motion and the motion of the suspended particles. For example, in a 

high fluid acceleration in the diverging region of a supersonic wind tunnel 

or in a rapid deceleration zone across a normal shock front, the velocity 

lag of the seed particles acting as the light scatterer would cause an 

appreciable error in the measurement of the fluid velocity employing a 

1 2 

Laser Doppler Velocimeter. ’ The size, the concentration, and the physical 
characteristics of the scattering aerosol are of paramount importance in the 
application of a Laser Doppler Velocimeter to fluid turbulence measurements. 
In general, the upper size limit is determined by particle inertia, and the 
lower size limit is determined by the increasing Brownian motion and the 
molecular slip in addition to the decreasing light scattering cross-section 
as particle size decreases. A polydisperse aerosol may cause an instrumental 
spread of the measured velocity distribution since a wide spectrum of aerosol 
size would result in a wide distribution of particle velocity lags. Thus, 
the scattering aerosol should be fairly monodisperse in size. The particle 
concentration should be sufficiently low so as to produce a negligible 
perturbation on the flow field, yet the concentration should provide the 
necessary amount of data in order that meaningful results can be obtained 
within the time of measurement. The physical properties of particles should 
be such that the aerosol is fairly stable, non-corrosive, non-toxic, and 
compatible with the physical properties of the fluid medium and the flow 
dynamics. Further, particle deposition on the optical windows should be 
minimum or should create minimum light scattering from the windows. The 
rate of generation and the method of introduction of the aerosol to the flow 
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medium should be controllable. 

The objectives of the present study are (1) to study different methods 
of generating aerosols having the desired physical properties, and (2) to 
characterize the behavior of aerosol particles in a flow field. Because of 
the nonlinear particle-fluid interactions involved, particularly when the 
particle Reynolds number is greater than one, emphasis is given to the 
experimental determination of the response of the particles in a gaseous flow 
field. An experimental study on the particle-fluid interaction may be 
performed by subjecting the seed particles to known sinusoidal fluid velocity 
fields or to known fluid velocity gradients. A possible method of generation 
of a known fluctuating or oscillatory flow field is the excitation of the 
fluid medium by a propagating acoustic field! The frequency response of a 
flow-borne particle can be measured by subjecting the particle to an acoustic 
field of known frequency. Since nonlinear couplings are involved, the 
particle response will be a function of the excitation frequency, particle 
Reynolds number, and particle relaxation time. The dynamic components of the 
acoustic accelerating force applied to the particle should correspond to 
the components of the fluid accelerating force that a particle will actually 
experience in the flow field to be investigated by using a Laser Doppler 
Velocimeter. 

This report discusses briefly the motion of aerosol particles in an 
acoustic field and some of the experimental arrangements for the particle 
response studies. It also describes two methods of generating a fairly mono- 
disperse scattering aerosol suitable for flow seeding for air turbulence 
measurements using Laser Doppler Velocimeters. Experimental data on small 
amplitude particle response are presented. 
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II. EQUATION OF MOTION OF A SMALL SUSPENDED PARTICLE 
A. Particle-Fluid Interaction 

For a small free spherical rigid particle in a locally uniform flow 
field, the force balance can be written from Newton's second law 

m p(dy dt ) s f d + F e (1) 

where is the mass of the particle, Vp is the particle velocity, F D is the 

fluid resistance force, and F g represents the summation of the forces acting 

on the particle other than the fluid resistance forces. Examples of the 

latter forces are: gravitational force, "buoyant" or "lift” force caused by 

pressure gradients as well as by the presence of a shear layer of fluid, 

3-5 

thermal forces and electrostatic forces. 

The fluid resistance force F^ can be expressed as 

F D= C D< 1/2 > r pV U g - V 1 U 9 • V (2) 

where C D is the drag coefficient, rp is the particle radius, p^ Is the 
average value of the fluid density, is the vectorial velocity of the 
fluid, and Vp Is the vectorial velocity of the particle. The drag 
coefficient, C Q , is a function of the particle Reynolds number, Re, which 
is given by 


RS “ d p V p9 5 9 /S <3) 

where v pg' lu g‘V’ 

and pg and n are the average values of density and viscosity of the fluid, 
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respectively, and d g is the particle diameter. 

For small Reynolds numbers, Re < 1 , where the inertial forces are small 
in comparison with the viscous forces in the medium, can be represented 


C D = 24/(ReC c ). (4) 

The above relationship indicates that the particle motion is in the Stokes' 
law regime. The term C represents the Cunningham correction factor given 
by 


C c * 1 + (2A*g/d p ) 


(5) 


where A = 1.25 + 0,44 exp( -1.09 d /2x ) 

r y 


( 6 ) 


5 

and x g is the mean free path of the gas molecules. For 700 <. Re <_ 2 X 10 , 
C D = 0.44 as was shown by Newton. For intermediate Reynolds numbers, 

1 2 

1 < Re < 700, several approximate empirical relationships are available 5 
which express C D in terms of Re, and it becomes difficult to give an exact 
analytical description of the flow field. 

If only one coordinate-component of U g and Vp are considered, and are 
denoted by u g and Vp, the equation of motion of a small spherical particle 
can be approximated in the following form:** 


f t 

(dVp/dt) + aVp = au g + b(du g /dt) + c l { (du g /dt ' - dv p /dt’)/(t - t') 1 ^ 2 }dt' 


where 


a = 36n/{(2pp + p g )d p 2 > 
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b « 3Pg/{2 Pp + p g ) 

and c * 18(p g n/nr)^ 2 /{(2pp + p g )dp) 


Tchen 6 solved this equation to obtain the ratio of the Lagrangian energy- 
spectrum function E (w) for the particle to the Lagrangian -energy spectrum 

r 

function for the fluid E g U). The relationship can be expressed as 

E p M/E g M - n + + f z Z M (8) 

where ' f.(„) = da _ 

(a + c/7tw/2) + (u + c/jTw/2 ) 


fgCw) = 


ui(a + c/tuij/ 2) (b - 1) 

(a + C/irw/2) + (o 3 + C/TTti)/2) 


<n is the angular frequency of turbulence components, and the constants a, b, 
and c are defined as above. Figure 1 shows the variation of E («)/£{ w) as 

r y 

a function of turbulence frequency for unit density spherical particles of 
1.0 and 0.3 micron-diameter calculated from Equation (8). 

When the particle Reynolds number Re < 1, and p p >> p g , the terms 
containing the constants b and c in Equation (7) can be neglected. Thus, if 
the particle motion is in the Stokes' law regime, Equation (8) can then be 
simplified as follows 


Ep(w)/E g U) = 1/(1 + Ap 2 ) (9) 

7 

Where x p is the dynamic relaxation time of the particle which can be 
expressed as 
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r p - (Va)C c = Pp d p 2 C c /18S (10) 

* 

The term C , the Cunningham correction factor, becomes significant when d 

C r 

is smaller than or comparable to A g , the molecular mean free path. The 
mean free path increases significantly as the flow approaches the free 

g 

molecular regime in high mach number flows. 

Equations (8) and (9) are based on linear fluid resistance law. In a 
low-Reynolds-number flow, the fluid resistance force is given by Stokes' law, 
i.e.. 


f d = 6tt ^ v p " U g ) (11) 

g 

and the heat transfer between the gas and a single particle is 

% - «V (T P - T g ) (12) 

where k is the thermal conductivity of the gas, T and T are the instantaneous 

I *7 

temperatures of the particle and of the surrounding fluid respectively. 
Equations (11) and (12) are valid for a spherical particle moving 5 with uniform 
velocity at a low Reynolds number in a viscous fluid. In most actual 
turbulent flow cases, temperature, density and dynamic viscosity of the fluid 
are variables and the fluid resistance force becomes nonlinear. 

Motion of suspended particles in an oscillating fluid can be studied by 
subjecting the aerosol to an acoustic plane wave. Since the energy of 
turbulence in a fluid flow is distributed over a spectrum of frequencies, 
experiments carried out at discrete frequencies can be related to the motion 
of particles in a turbulent fluid if the fluid turbulence conditions are 
simulated. Since nonlinear fluid resistance forces are involved, both 
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amplitude and frequency of the fluctuating fluid velocity field under test 
conditions should be comparable to the turbulent velocity components over 
the spectrum of frequencies. 

Figure 2 illustrates typical values of root mean square turbulence 
velocity components as a function of frequency of turbulence in a subsonic 
wind tunnel. The mean velocity is approximately 300 meters per second. 
Curve No. 1, shown as an unbroken line, represents turbulent velocity components 
in a wind tunnel of throat area approximately 30 cm X 30 cm. Curve No. 2, 
shown as an unbroken line, represents similar quantities for isotropic 
turbulence in a 25 cm-diameter wind tunnel. The curves shown in broken 
lines represent corresponding rms Reynolds numbers of a spherical, unit 
density, one micron-diameter particle suspended in the fluid medium. The 
fidelity of seeded particles in tracking the fluid velocity fluctuations in 
such a flow field can be experimentally determined by measuring the ratio of 
the particle velocity amplitude to the fluid velocity amplitude while the 
particles are in an oscillatory flow field of comparable velocity amplitude 
at the frequency of excitation, 

B. Tracking Fidelity of Aerosols in an Oscillatory Flow Field 

The equation of motion of a particle of radius rp in an oscillatory 

7 

flow field can be written in the form 

m p( dv p /dt ) * ~ 6 ^ r p( v p ■ V 0 3 ) 

+ 6Trnr p (wr p 2 /2v) 1/f2 [(u g - v p ) + (l/w)d(u g - v p )/dt 

+ <2/3w)Ur p 2 /2v) 1/2 (du g /dt)j 

where v is the kinematic viscosity of the gas and u is the angular frequency 
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of oscillation. The second term of Equation (13) represents deviation from 
Stokes' law regime. If the nonlinear terms in Equations (7) and (13) are to 
be made negligible, the following constraints must be satisfied: 


WT d- 1 
Ur p 2 /2v) 1/2 « 1 

(p g /p p ) « 1 

Og <. 1.5 ( 14 ) 

Re << 1 
K n < 0.2 

nipO/pg < 0.01 


Under these conditions of linear particle-fluid interactions, a high tracking 
fidelity of a flow borne aerosol in a subsonic flow field can be assumed. 

The above conditions are not mutually independent but do set criteria 
relating to the physically meaningful quantities. The first three criteria 
imply that the flow Is in Stokes' regime. These conditions also imply that 
the dynamic relaxation time T p and thermal relaxation time t.j. of the 
particle are small. The latter quantity is given by 


* p ,V r p 2e p /3 '’ c P 


(15) 


where P r is the Prandtl number which is equal to Cpfl/k, c p ' is the specific 
heat of the particle, c p is the specific heat of gas at constant pressure 
and k is the thermal conductivity of the gas. The fourth condition implies 
that the seeded aerosol should be fairly monodisperse and the standard 
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deviation of the particle size spectrum is close to one. Because of 
inherent difficulties in generating a true monodisperse aerosol and because 
of the inherent unstable nature of aerosols, a maximum value for the standard 
deviation (a ) of 1.5 is chosen as a practical limit. A larger value of a 
may cause appreciable error in turbulence measurement because of the possible 
wide distribution of particle velocity lags. The fifth constraint is 
satisfied in most actual subsonic flow cases where r p < 2 microns, 

P p * 1 gram/cc and the frequency of turbulence does not exceed 100 kHz. 

Figures 1 and 2 indicate that for d p = 1 micron, the gas-particle relative 
velocities and rms Reynolds numbers remain small in subsonic turbulence. 

The requirement that the Knudsen number (K R = x^/dp) be smaller than 0.2 
implies that for air at STP, dp > 0.3 microns, so that the molecular slip is 
small. The last constraint indicates that particulate mass loading ratio is 
small; thus the perturbation of the flow field will be negligible. The lower 
limit of particulate mass loading is set by the Doppler signal dropouts. 

In discussing the above constraints, it was assumed that the scattering 
particles are uncharged, inert, and are not volatile in the surrounding fluid. 

Figure 3 shows the calculated maximum values of the frequency of the 
oscillating flow field and of the diameter of unit density spherical 
scatterers for v p /u g = 0.707 (m d = 1.0) and for v p /u g = 0,980 (d p 4>" = 0.02). 
Assuming that a tracking fidelity of 98 percent is desired in a turbulence 
measurement where the spectral components do not exceed 100 kHz, d p AT = 0.02, 

r 

and f = 10 Hz, the value of the particle diameter, d p , is 0.25 microns for 
98 percent fidelity. Considering the constraints in Equation (14) that 
K n < 0.2, i.e., d p > 0.3 microns for air at STP, a possible optimum choice 
of tracer particle diameter would probably be 0.3 microns with p p » l.o for 
air turbulence measurements. 
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C. Motion of a Particle in an Acoustic Field 

When the constraints given in Equation (14) are satisfied, the motion 
of an oscillating particle in an acoustic field, can be written in the 
form 


T n(dv,Vdt) + v o = u o sin wt (16) 

p' p p g 

Equation (16) represents one dimensional motion of a particle in an acoustic 
field of angular frequency o». Equation (16) has a solution in the form 


» (u g sinUt + <j>)/( 1 + o> 2 x p 2 )^ 2 } 
+ (Tp U g w/O + W 2 T p 2 )}e p 


where the phase shift $ is given by 


» tan'^uT. 


07 ) 


(18) 


The second term of Equation (17) represents the transient response of the 
particle. When t » t , we can write the steady state solution for the 

r 

ratio of the amplitude of the particle velocity to the amplitude of the 
fluid velocity as 

(v p /u g ) - 1/(1 + <A p 2 ) 1/2 (19) 

The ratio (Vp/u g ) is often referred to as the degree of entrainment and is 
denoted by Equation (19) indicates that the fidelity with which an 
aerosol particle follows the fluid motion increases with decreasing particle 
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diameter, density, and excitation frequency and with increasing viscosity of 
the medium. 

The velocity amplitude of an element of fluid subjected to an acoustic 

2 

plane wave of intensity I in watts/meter , is given by 

“g - (2I/p gCg ) 1/2 (20) 

where c is the velocity of the propagating acoustic wave. The amplitude 
of displacement is given by 

A g = u g /» (21) 

where u is the angular frequency of excitation. Since acoustic intensity 

-12 2 

levels are generally measured in terms of dB referred to 10 watts/meter , 
the intensity level IL is given by 

IL - 10 log,. 1 l a . "At t s/m 2 , mtu/ J (22) 

- ,u I Q in watts/nr 0 

where I - 10" watts/meter . Figure 4 shows calculated relationships 
between acoustic intensity levels in dB and fluid velocity amplitude u^ 

(solid line). The velocity amplitude, v p , of a suspended particle can be 
calculated from Equation (19). Using Equations (3) and (18) the particle 
Reynolds numbers were calculated for d p = 1, 5 and 10 microns; p p * 1.0 
gram/cc; and f = 10 and 100 kHz. The calculated values of Re are plotted 
in Figure 4 (broken lines). Air at 760 mm Hg at 20°C is assumed as the 
fluid medium. 

12 13 

Several experimental verifications of Equation (19) have been attempted * 
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most of which employed determination of the amplitude of particle oscillation 
using photographic means. Although there was general agreement between the 
experimental data and the expected particle response, the experimental results 
scatter widely and the validity of the analysis was more qualitative than 
quantitative. A Laser Doppler Velocimeter provides an excellent method of 
measurement of suspended particle response in acoustic fields. 

In Figure 4, a linear coupling of acoustic excitation with the medium 
is assumed. In practice, the waveform of the propagating wave becomes 
distorted at a distance from the source, particularly when the acoustic 
intensity is high. If nonlinear terns in particle-fluid interaction are 
present, particle motion may become further distorted. However, the 
approximate values of the particle Reynolds numbers plotted in Figure 4 
indicate that the acoustic intensity level required to subject an aerosol 
(with d = 1 micron, p = 1.0 gram/cc) to a value of Re _> 1 would be very 

r r 

high. Several particle difficulties limit the generation of such a highly 
intense acoustic field in gases. However, an oscillatory flow field of 
high velocity amplitude can be established in a resonant air column such as 
a siren tuned cavity. Gal ton whistle, etc. 

Particle-fluid interaction with a particle Reynolds number Re > 1, may 
be realized in the following flow conditions: (1) particles suspended in a 

fluid acceleration field passing a throat, (2) flow of particles through a 
standing shockwave, where particles experience a large deceleration field 
and (3) propagation of a shockwave through a stationary fluid medium 
containing particles. 

III. GENERATION OF A MONODISPERSE AEROSOL 


Two of the most desirable features of seed particles in Laser Doppler 
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Velocimeters are that (1) the particles must be an efficient scatterer at 
the wavelength of laser radiation and (2) flow tracking fidelity must be 
high. As discussed in the earlier section, a monod is perse submicron aerosol 
of solid particles or of liquid droplets seems to satisfy these requirements. 
Liquid droplets of very low vapor pressure have an advantage over the solid 
particles in that their deposition on the optical windows results in a 
formation of a liquid film that does not cause a severe light scattering 
problem from the window surfaces. Since aerosol coagulates spontaneously, 
the coagulated liquid droplets remain spherical. However, solid particles 
with a high melting point must be used where temperatures of the fluid 
medium are high. 

A. Aerosol of Submicron Liquid Droplets 

Laskin-type atomizers^ can be used to generate DOP (p = 0,98; refractive 
index, m = 1.4 at X = 5000A; vapor pressure = 1.48 X 10 ^ mm Hg at 25°C) 
aerosol as well as Dow Corning 200 Silicone fluid (p = 0.96; m = 1.4 at 
x = 5000A; viscosity 50 and 100 centistokes) aerosols. As a rule, pneumatic 
atomization produces a polydisperse mist. A typical cumulative size distribu- 
tion of DOP aerosol 15 produced by a Laskin atmoizer is shown in Figure 5. 
Figure 5A shows a typical photomicrograph of a sample of the aerosol generated 
by the Laskin atmoizer. 

Generation of a monodisperse condensation aerosol using a heat-exchanger 
type generator 15 is advantageous compared to the Laskin atomizer. A 
Rapaport-Weinstock generator was used for generating fairly monodisperse 
DOP aerosols. Aerosols with mean diameter varying from 0.1 to 1.5 microns 
could be produced. A schematic diagram of the generator is shown in 
Figure 6. A dilute solution of DOP in alcohol was atomized and the poly- 
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disperse mist was evaporated in the heated column. The vapor finally 
condensed on the impurity nuclei to form a monodisperse aerosol. By using 
different nozzle sizes and by controlling the dilution ratio of DOP to 
alcohol, the droplet size could be varied from 0.1 to 1.5 microns in diameter. 
Two typical photomicrographs of DOP aerosol are shown in Figures 7 and 7b. The 
droplets were electrostatically deposited on a microscope slide surface 
deactivated with 3M L-1668 chemicals. The standard deviation, o g , of the 
particle size distribution can be maintained below 1.5 

B. Aerosols of Monodisperse Solid Particles 

Monodisperse aerosols containing solid spherical particles were 
generated by atomizing an alcohol suspension of polystyrene latex spheres 
(p p = 1.05, m = 1.58 at \ = 5000A). Uniform (one standard deviation within 
0.005 microns) polystyrene latex particles in the size range 0.1 to 1.0 
micron were used. These particles were obtained from Dow Diagnostics. The 
aerosol generator is shown in Figure 8. A lucite tube, 1/4" I.D. and 
1/2" O.D., with 6 radial holes of 0.0135" diameter was used as an atomizer. 

The bottom end of the tube was closed. The open end was connected to a 
pressurized Ng tank through a pressure regulator and an absolute filter. 

The radial holes were positioned so that the nitrogen gas jets skim the 
surface of the suspension. A dilute suspension was used to produce an 
aerosol in which most of the. particles were singlets (latex particle). It 
was found that the latex particles in an alcohol suspension have less 
tendency to coagulate when the solution concentration is altered. Also, 
alcohol film on the particles evaporates faster than the water film. Thus 
an alcohol suspension was used in the latex particle aerosol generator. 

Figures 9 and 9a show two photomicrographs of the generated aerosol of 0.500 
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1.011 ml cron -diameter particles. The ratio of single particles to all other 
particles, such as doublets, triplets, etc. was found to be high. 


IV. EXPERIMENTAL RESULTS 
A. Experimental Setup 

An experimental arrangement for the measurement of the frequency response 
of a scattering aerosol is shown in Figure 10. A 500 mW argon-ion laser beam 
of wavelength 4880A is incident on the outermost track of an encoder disc 
having a spatial frequency of 100 lines/mm. The disc is rotated at a - 
sufficient speed by a synchronous motor to permit a frequency translation* 6 
of 2 MHz of the Doppler signal. The test aerosol is suspended in an acoustic 
excitation chamber. Air at room temperature was used as the fluid medium. 

The fluid medium containing the aerosol was excited by acoustic horns. The 
frequency of excitation could be varied at discrete steps from 100 Hz to 
100 kHz. The Doppler signal was demodulated by a phase-locked loop 
demodulator using a Signetics NE 561B PLL circuit. The acoustic pressure 
levels at the point of velocity measurements were measured by using a 
calibrated Bruel and Kjaer condenser microphone (1/4 inch diameter. Type 
4135) followed by a preamplifier (B&K model 2618). The demodulated Doppler 
signal and the microphone output were displayed on a dual beam oscilloscope. 
Velocity amplitudes of u^ and Vp were calculated from the oscilloscope 
displ ay. 

Acoustic excitation in the frequency range below 15 kHz was carried 
out using loud speakers and tweeters. Sound pressure levels in the range 
of 100 to 120 dB were used at different resonant frequencies of the chamber. 
For the transmission of ultrasonic vibrations into the fluid medium, flexural 
vibrating plates were used. A typical transducer driven by a piezoelectric 
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crystal is shown in Figure 11. The transducer consists of a circular, 
flexural ly vibrating, free-edge plate with stepped thicknesses. The plate 
is driven at its center and it can vibrate at several resonant frequencies 
corresponding to different flexural modes of vibration. A number of 
ultrasonic radiators were made with different resonant frequencies. The 

t 

, measured output acoustic levels in air at a distance of 3 cm from one of 
the ultrasonic radiators is shown in Figure 12. 

During the acoustic excitation, standing waves in both longitudinal 
and transverse modes were formed. Since only a one-dimensional particle 
oscillation was measured with the Laser Doppler Velocimeter, a correction 
factor was needed for obtaining an exact coincidence of v / u = 1.0 for 

r y 

dp = 0.176 microns at f = 20 kHz as is expected from Equation (19). 

B. Frequency Response of Aerosol 

The method of determining the particle response in an oscillatory fluid 
velocity field was as follows: The test aerosol was introduced inside the 

acoustic excitation chamber by connecting the input of the chamber to the 
aerosol generator while the other part of the chamber was connected to a 
vacuum line. After a suitable concentration of aerosol was obtained inside 
the chamber, both the ports were closed. The acoustic exciter was then 
energized and the acoustic pressure level at a point close to the laser beam 
crossing was measured by the microphone. The velocity oscillations of the 
aerosol were measured by the Laser Doppler Velocimeter. Both the modulating 
signal (microphone output) and the demodulated signal (Doppler signal 
processor output) were displayed on a dual beam oscilloscope. A sample of 
aerosol was taken from the chamber and its photomicrograph was analyzed to 
determine the aerosol size distribution. The process was repeated for each 
test aerosol. 
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Figures 13 through 19 show some typical oscilloscope displays of 
modulating and demodulated signals for different particle size and for 
different frequencies of acoustic excitation. Figures 13 and 14 show the 
response of a polydisperse DOP aerosol in acoustic fields. The aerosol was 
generated by a Laskin atomizer. As the modulating frequency increased from 
5 kHz to 50 kHz, the velocity amplitudes of the aerosol showed a velocity 
distribution as is evident from the spread in amplitude of the displayed 
signal. Such a velocity spread is expected at high frequencies since the 
larger particles lag the smaller ones. Figure 15 shows a similar display 
of the reponse of a polydisperse DOP aerosol using DISA Type 55L35 Frequency 
Tracker as the Doppler signal demodulator. The tracker was obtained on 
loan from DISA Electronics, New Jersey, in order to determine whether or 
not the signal-to-noise ratio of the demodulated output can be further 
improved by using a frequency feedback system such as is used in the tracker. 
No improvement in SNR was observed compared to the phase-locked-loop system 
in the present application. A bandpass filter having a narrow bandwidth 
centered around the modulating frequency was found necessary for Doppler 
signal demodulation. In the present work, a phase-locked-loop FM demodulator 
and tunable bandpass filter were used for Doppler signal processing. 

Figure 16 shows the response of a fairly monodisperse DOP aerosol at 
the modulation frequency of 20.84 and 51,46 kHz. The aerosol was generated 
by using the Rapaport-Heinstock aerosol generator. A comparison of Figure 
16 with Figures 13 and 14 shows that the distribution of velocity lag is much 
narrower in the case of the DOP aerosol generated by the Rapaport-Weinstock 
aerosol generator. The narrower size distribution of the aerosol is also 
evident from the photomicrograph shown in Figure 7A when the same is compared 
with the size distribution shown in Figure 5A. 
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Figures 17 through 19 show the response of monodisperse aerosols in 
acoustic fields. The distribution of velocity lag for such test aerosols 
was found negligible which also indicates that the aerosol particles are 
mostly singlets and the coagulation rate is small. 

Acoustic excitation levels at the point of velocity measurement were 
maintained at a constant level for a given frequency of excitation. The 
particle velocity amplitudes for different size particles can thus be 
directly compared from the oscilloscope display of the demodulated signal. 

Thus the ratio of the velocity amplitudes was not affected by the attenuation 
of the acoustic wave caused by the presence of the aerosol. An appreciable 
attenuation was found when the mass loading ratio (m n/p ) was high and the 
product ojtp was close to one. 

Table 1 shows the ratio of the measured particle velocity amplitude for 
the test aerosol to the measured particle velocity amplitude for particles 
of 0.176 micron-diameter (p = 1.05). The choice of 0.176 micron-diameter 

r 

particles as a reference for u^ was made since these particles are expected 
to follow an oscillatory flow field without any appreciable velocity lag in 
the frequency range DC to 100 kHz. For example v p / u g ~ 0.994 for f = 84.0 kHz 
(from Equation 19). The relative amplitude ratio (v p /v Q 175) was preferred 
since the measurement of u g from the microphone output may have the following 
sources of error: (1) the microphone is not directional' and thus it cannot 

measure the one dimensional component u g alone if other orthogonal components 
are present, and (2) the size of the microphone becomes comparable to the 
wavelength of acoustic radiation at the higher frequencies. Figure 20 shows 
a plot of the theoretically calculated ratio of relative particle velocity 
amplitudes. Experimental data points show a close fit to the expected 
values. 

Figure 21 shows the entrainment ratio (v /u ) as a function of r w 1 / 2 

1 V I 



- 20 - 


(a factor related to the square root of the nondimensional frequency wT p ). 

In this figure, the experimental data were taken for d p = 0.176, 0.500, 

0.760 and 1,011 microns polystyrene particles and for f = 8.56, 37.80, 58.00 
and 02.80 kHz. In this experiment a different horn was used for acoustic 
excitation. Although the 0.176 micron-diameter particles follow the 
velocity fluctuations with a high fidelity, the optimum choice would probably 
be 0.30 micron-diameter tracer particles because of the larger light 
scattering crosssection of 0.3 micron particles. 

The average acoustic intensity level at the point of velocity measure- 
ment was approximately 120 dB which corresponds to a value of u^ = 5.5 cm/sec. 
Thus the experimental data presented here are strictly valid for low amplitude 
oscillation of particle velocity. However, if the velocity amplitude is 
increased further to simulate the expected turbulence component as shown in 
Figure 2, the constraints of Equation (14) would still be valid. As a result 
no appreciable departure from the results shown in Figures 20 and 21 is 
expected. An experimental verification of particle fidelity under simulated 
conditions close to the actual turbulent flow field is now underway. 

C, LDV Signal Demodulation 

During the present experiments, it was observed that laser Doppler 

signal processing at a modulating frequency above a few kHz becomes a 

difficult problem. Doppler signal demodulators, used in the present work, 

show that when the frequency (f) of a fluid velocity oscillation is 

increased above a certain value, the output SNR of the signal degrades to a 

point where no meaningful data can be taken. The degradation of output SNR 

is probably caused by the fact that as f is increased the modulation index 
1 2 

g(Af/f) decreases and the output SNR (SNR aB ) decreases. When the SNR is 
below the threshold value of FM demodulation, the performance of the analog 
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tracking systems becomes very poor. The degration is further accentuated 
by signal dropouts. SNR can be improved significantly if the frequency of 
oscillation is known. In such a case a narrow bandpass filter can be used 
to filter the noise outside the pass band. A tunable bandpass filter was 
used in the present set of experiments. In actual turbulence analysis, the 
processor must have a demodulated signal bandwidth DC to 100 kHz. A digital 
counter type or a photon correlator type of Doppler signal demodulator may 
be suitable in such applications. 


V. CONCLUSIONS 

An analysis of the particle-fluid interaction and scattering properties 
of the aerosol indicates that a fairly monodisperse aerosol of 0.3 micron- 
diameter is an optimum choice for tracer particles in subsonic turbulence 
studies employing a Laser Doppler Velocimeter, Experimental results presented 
here seem to confirm this analysis. Methods of generating a submicron aerosol 
for flow seeding of inert gases at ambient temperature are discussed and 
experimental data are presented. 
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TABLE 1 


Correlation of Theoretically Calculated and 
Experimentally Observed Velocity Amplitude Ratios 
(Note: Vq ^g/u ~ Q»994 0 84 kHz) 


v p /v 0.176 


Frequency 

of 

Oscillation 

kHz 

Particle 
Diameter 
in Microns 

Experimental 

Theoretically 
Calculated 
(From Stokes Law) 

4.99 

0.500 

1.00 

0.99 

4.99 

0.760 

1.00 

0.99 

4.99 

1.011 

1.00 

0.99 

4.99 

2.020 

0.88 

0.91 

20.84 

0.500 

1.00 

0.99 

20.84 

0.760 

0.97 

0.96 

20.84 

1.011 

0.92 

0.89 

20.84 

2.020 

0.43 

0.48 

32.15 

0.500 

1.00 

0.98 

32.15 

0.760 

0.90 

0.92 

32.15 

1.011 

0.79 

0.80 

32.15 

2.020 

0.32 

0.33 

51.46 

0.500 

0.90 

0.95 

51.46 

0.760 

0.73 

0.82 

51.46 

1.011 

0.60 

0.64 

51.46 

2.020 

0.19 

0.21 

84.80 

0.500 

0.82 

0.88 

84.80 

0.760 

0.66 

0.65 

84.80 

1.011 

0.45 

0.45 

84.80 

2.020 

0.13 

0.13 




FIGURE 1 

DECREASE OF ENERGY SPECTRUM RATIO E ( U )/E (w) WITH INCREASING 

p g 

TURBULENCE FREQUENCY AS A FUNCTION OF PARTICLE SIZE 
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FREQUENCY OF TURBULENCE *- 

FIGURE 2 

TYPICAL RMS TURBULENT VELOCITY COMPONENTS OF SUBSONIC FLUID FLOW IN A WIND TUNNEL (U 


300 m/sec) AND 


THE CALCULATED REYNOLDS NUMBER OF A UNIT DENSITY SPHERICAL PARTICLE PLOTTED AS A FUNCTION OF TURBULENCE FREQUENCY 


TURBULENCE FREQUENCY (f) 







FIGURE 4 


CALCULATED VARIATION OF THE FLUID VELOCITY AMPLITUDE (u ) AND 

THE ASSOCIATED PARTICLE REYNOLDS NUMBER IN AN ACOUSTIC FIELD 
PLOTTED AS A FUNCTION OF ACOUSTIC INTENSITY LEVEL 
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PERCENT OF MASS ASSOCIATED WITH PARTICLES OF DIAMETER LESS 

THAN THAT INDICATED 
FIGURE 5 


CUMULATIVE TYPICAL SIZE DISTRIBUTION CURVE FOR DOP 
AEROSOL GENERATED BY LASKIN ATOMIZER 





FIGURE 5A 

A PHOTOMICROGRAPH OF DOP AEROSOL GENERATED BY A 
LASKIN ATOMIZER (LINEAR MAGNIFICATION = 2090X) 


DOP AEROSOL 



no v 

60 Hz 
SUPPLY 


RAPAPORT- WE INSTOCK AEROSOL GENERATOR 





‘p'v, 


o 

CO 


o 

Oo o 

D©° o° 

O 


O 

O 0 o O 
o o 
o 


o o 


Oo W o 
°o O 0 

o <P° 


Oo° 


c 

c 


0°0 


o 

o 

O c Co 



A r\ 


° o ° 

o o o 

o o ° 

o 

o o 

O ° o O 


o 

O o 

C Prs C 


c 

O o° 

'o°° 

o 


c 

c o 

o 


o 

c 


Qc 


FIGURE 7 

A PHOTOMICROGRAPH OF OOP AEROSOL GENERATED BY RAPAPORT-WE INSTOCK 
GENERATOR USING A 25 PERCENT CONCENTRATION OF DOP (LINEAR MAGNIFICATION = 2090X) 
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FIGURE 7A 

A PHOTOMICROGRAPH OF UOP AEROSOL GENERATED BY RAPAPORT-WEINSTOCK 
GENERATOR USING A 5 PERCENT CONCENTRATION OF DOP (LINEAR MAGNIFICATION = 2090X) 


FILTERED 



FIGURE 8 


POLYSTYRENE LATEX PARTICLE AEROSOL GENERATOR 



FIGURE 9 

A PHOTOMICROGRAPH OF 0.500 MICRON-DIAMETER POLYSTYRENE 
LATEX AEROSOL (LINEAR MAGNIFICATION = 2090X) 
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FIGURE 9A 

A PHOTOMICROGRAPH OF 1.011 MICRON-DIAMETER POLYSTYRENE 
LATEX AEROSOL (LINEAR MAGNIFICATION = 2090X) 
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FIGURE 10 I 

AN EXPERIMENTAL ARRANGEMENT FOR MEASURING THE FREQUENCY RESPONSE OF 
SCATTERING AEROSOL OSCILLATING IN AN ACOUSTIC FIELD 








FIGURE 11 

SCHEMATIC DIAGRAM OF A CIRCULAR PLATE RADIATOR 
(A = WAVELENGTH OF SOUND IN THE RADIATOR MATERIAL) 




INXfNSITY LEVEL (dB) 



FIGURE 12 

ULTRASONIC INTENSITY LEVELS MEASURED AT A DISTANCE OF 3 CENTIMETERS FROM 
THE TOP OF THE RADIATING PLATE AS A FUNCTION OF EXCITATION FREQUENCY. 
INPUT VOLTAGE 600 VOLTS PEAK-TO-PEAK, RESONANCE FREQUENCY = 37.8kHZ. 


(a) Frequency of Excitation: 4.99 kHz 
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(b) Frequency of Excitation 20.84 kHz 
FIGURE 13 

RESPONSE OF A POLYD1SPERSE OOP AEROSOL ( ffg > 2.0, AVERAGE DIAMETER 
1.03 MICRONS, CONCENTRATION = 10 7 PARTICLES/cc) IN ACOUSTIC FIELDS 





Aerosol 

Response 


Microphone 

Output 



a. Frequency of Excitation: 32*15 kHz 
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b. Frequency of Excitation: 51.46 kHz 
FIGURE 14 

RESPONSE OF A PQLYDISPERSE OOP AEROSOL (o g >2.0, AVERAGE DIAMETER 

1.03 MICRONS, CONCENTRATION = 10 7 PARTICLES/cc) IN ACOUSTIC FIELDS. 
PARTICLE SIZE RANGE: 0.3 - 2.0 MICRONS 
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a. Frequency of Excitation: 20.84 kHz 
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b. Frequency of Excitation: 32.15 kHz 
FIGURE 15 


POLYDISPERSE DOP AEROSOL RESPONSE IN ACOUSTIC FIELDS MEASURED 
BY A DISA 55L20 DOPPLER SIGNAL PROCESSOR 
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b. Frequency of Excitation: 51.46 kHz 
FIGURE 16 

FAIRLY MONODISPERSE DOP AEROSOL (o g = 1.4, AVERAGE DIAMETER » 0.6 MICRONS, 
CONCENTRATION = 10 6 PARTICLES/ cc) RESPONSE IN ACOUSTIC FIELDS 
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b. Frequency of Excitation: 51.46 kWl 

FIGURE 17 

MONOOISPERSE POLYSTYRENE AEROSOL (o g = 1.0, PARTICLE OJAKETER » 1.01 MCRONS, 
CONCENTRATION = I0 5 PARTICLES/ccJ RESPONSE IN ACOUSTIC FIELDS 
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a. Frequency of Excitation: 20.84 kHz 
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b. Frequency of Excitation: 51.46 kHz 

FIGURE 18 

MQN0D1SPERSE POLYSTYRENE AEROSOL (<r g = 1.0, PARTICLE DIAMETER - 0.50 MICRONS, 
CONCENTRATION - 10^ PARTICLES/cc) RESPONSE IN ACOUSTIC FIELOS 




Aerosol 

Response 


Microphone 

Output 


EHSHMi 

wamBgmaa 

iiMiiin 

wMnmmwmwm 
nNinifiinm 

mmmmm 

a. Frequency of Excitation: 20.84 kHz 


Aerosol 

Response 
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Output 



b. Frequency of Excitation: 51.46 kHz 

FIGURE 19 

M0N0D IS PERSE POLYSTYRENE AEROSOL (a g = 1.0, PARTICLE DIAMETER = 0.17 MICRONS, 
CONCENTRATION = 10 5 PARTICLES/cc) RESPONSE IN ACOUSTIC FIELDS 
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FIGURE 20 

COMPARISON OF EXPERIMENTALLY MEASURED AND THEORETICALLY CALCULATED VELOCITY 
AMPLITUDE RATIOS FOR MONODISPERSE PARTICLES IN THE SIZE RANGE OF 
0.176 TO 2.02 MICRON-DIAMETER AND FOR THE FREQUENCY OF 
OSCILLATION IN THE RANGE 5 TO 85 kHZ 
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(a) Frequency of Excitation: 4.99 kHz 
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(b) Frequency of Excitation 20.84 kHz 
FIGURE 13 

RESPONSE OF A P0LY0ISPERSE OOP AEROSOL (<r g > 2.0, AVERAGE DIAMETER 
1.03 MICRONS, CONCENTRATION = 10 7 PARTICLES/cc) IN ACOUSTIC FIELDS 
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a. Frequency of Excitation: 32.15 kHz 
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b. Frequency of Excitation: 51.46 kHz 


FIGURE 14 

RESPONSE OF A POLYDISPERSE OOP AEROSOL (o g > 2.0, AVERAGE DIAMETER 

1.03 MICRONS, CONCENTRATION = 10 7 PARTICLES/cc) IN ACOUSTIC FIELDS. 
PARTICLE SIZE RANGE: 0.3 - 2.0 MICRONS 
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a. Frequency of Excitation: 20.84 kHz 



t. Frequency of Excitation: 32.15 kHz 
FIGURE 15 

POLYDISPERSE OOP AEROSOL RESPONSE IN ACOUSTIC FIELDS MEASURED 
BY A DISA 55L20 DOPPLER SIGNAL PROCESSOR 
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b. Frequency of Excitation: 51.46 kHz 
FIGURE 16 

FAIRLY IIONODISPERSE DOP AEROSOL (o g = 1.4, AVERAGE DIAMETER = 0.6 MICRONS, 
CONCENTRATION = 10 6 PARTICLES/ cc) RESPONSE IN ACOUSTIC FIELDS 
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a. Frequency of Excitation: 20.84 It Hz 



Aerosol 

Response 


Microphone 

Output 



L. Frequency of Excitation: 51.46 kHz 

FIGURE 17 

MONOD IS PERSE POLYSTYRENE AEROSOL (o g = 1.0, PARTICLE DIAMETER « 1.01 MICRONS, 
CONCENTRATION = IO 5 PARTICLES/cc) RESPONSE IN ACOUSTIC FIELDS 
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a. Frequency of Excitation: 20.84 kHz 
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b. Frequency of Excitation: 51.46 kHz 

FIGURE 18 

MONODISPERSt POLYSTYRENE AEROSOL (a g = 1.0, PARTICLE DIAMETER « 0.50 MICRONS, 
CONCENTRATION = 10 5 PARTICLES/cc) RESPONSE IN ACOUSTIC FIELDS 
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a. Frequency of Excitation: 20.84 kHz 
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b. Frequency of Excitation: 51.46 kHz 

FIGURE 19 

MONODISPERSE POLYSTYRENE ALROSOL ( a ^ 1.0, PARTICLE DIAMETER = 0.17 MICRONS, 
CONCENTRATION = 10 5 PARTICLES/cc) RESPONSE IN ACOUSTIC FIELDS 


